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Ruan 等 四 研究 者 利用 CHAMP (CHAllenging Min- 
isatellite Payload) 资 料 详细 研 究 了 MDM 特 征 ， 指 
出 它 与 太阳 活动 和 季节 有 关 , 并 用 物理 模型 研究 了 
MDM 的 形成 机 制 四. 
自 2001 年 以 来 , 3 颗 极 轨 卫 星 CHAMP GRACE- 
A (Gravity Recovery and Climate Experiment- 
A) 和 SWARM-C (The Earth's Magnetic Field and 
Environment Explorers-C) 以 较 高 的 时 空 分 辩 率 测 
量 高 层 大 气 密度 . 由 于 极 轨道 的 特点 ， 资 料 履 盖 
全 球 纬度 , 同时 在 纬度 70° 以 内 同一 上 升 段 (或 下 降 
段 ) 几 乎 对 应 同一 地 方 时 ， 上 升 和 下 降 段 分 别处 于 
日 照 区 和 阴影 区 ， 且 相隔 12 h( 有 时 恰好 处 于 晨 
线 轨道 面 ). 当地 方 时 接近 时 , 对 不 同 纬度 的 连续 测 
量 有 利于 研究 大 气 密 度 在 夜间 中 低 纬度 的 MDM 现 
象 . 相 比 于 过 去 的 研究 , 本 文选 择 以 上 3 颗 卫 星 的 资 
料 , 时 间 跨 度 更 长 (2001--2021 年 )， 通 过 多 星 对 比 
和 大 样本 统计 研究 , 可 以 更 清晰 地 证 认 MDM 特 征 ， 
获得 具有 统计 意义 的 定量 结论 . 

文章 首先 给 出 提取 MDM 特 征 的 方法 ,并 分 析 
地 方 时 、 高 度 、 辐 射 水 平 、 季 节 等 因素 对 MDM 
特征 的 影响 ,再 用 DTM2000 (Drag Temperature 
Model 2000)lsl, NRLMSISE00 (US Naval Research 
Laboratory, Mass Spectrometer and Incoherent 
Scatter radar)? fllJB2008 (Jacchia-Bowman 2008 
model)! 3; 339i E. Æ fir fe ie BE MD M3SE 47 BEA, 
并 与 观测 对 比 研究 . 最 后 , 尝试 采用 更 精细 的 经 验 


ftp://swarm-diss.eo.esa.int), 时 间 分 辨 率 为 10s. 由 
于 是 极 轨 卫星 , 每 一 圈 的 资料 都 能 覆盖 士 85" 以 内 的 
地 理 纬度 , 约 4 个 月 完整 遍历 一 次 0-360? 经 度 . 时 间 
范围 方面 ， 选 用 CHAMP 卫 星 2001 年 5 月 至 2008 年 
12 月 、GRACE-A 卫 星 2002 年 4 月 至 2014 年 12 月 、 
SWARM-C 卫 星 2014 年 12 月 至 2021 年 5 月 的 观测 资 
时 间 跨 度 近 20 yr. 地 磁 爆 发 会 在 高 纬度 地 区 产 
m 沿 着 纬度 方向 传播 上 ,造成 不 同 
密度 的 短期 涨 落 . 因此 , 为 了 研究 MDM 现 象 ， 
需 选 择 地 磁 平 静 期 资料 进行 研究 , 判 据 是 : 某 一 天 
平均 地 磁 指 数 Ap 平 均值 小 于 10, 且 前 一 天 3 h 地 磁 
指数 Ap 小 于 等 于 15. 经 筛选 , CHAMP、GRACE-A 
和 SWARM-C 分 别 有 913 d、1873 d、1078 d 的 资 
BL, 相互 之 间 有 时 间 重 登 ,， 有 利于 做 对 比 研 究 . 此 
外 ， 为 避免 高 度 起 伏 对 密度 的 影响 , 用 文献 [12] 给 
出 的 方法 , 将 每 颗 星 测量 的 大 气 密度 按 线性 比例 规 
律 归 算 至 其 所 在 的 平均 高 度 : CHAMP 为 360 km、 
GRACE-A 为 480 km、SWARM-C 为 450 km. 


2.2 ” MDM 特征 提取 

首先 以 天 分 组 , 挑选 轨道 落 在 阴影 区 的 大 气 密 
FE, 计算 南北 纬 70° 之 间 的 大 气 密 度 平均 值 作为 基 
准 密度 ; 然后 以 10° 为 宽度 ,计算 不 同 纬度 区 间 的 
密度 平均 值 , 再 与 基准 密度 求 比值 (下 文 简称 “密度 
EE"), 可 观察 每 一 天 大 气 密 度 随 纬度 的 连续 起 伏 . 
为 获得 统计 特征 ,将 多 天 资料 按照 不 同 策略 分 组 


建 模 函 数 , 对 MDM 观 测 特征 进行 拟 合 , 为 大 气 模型 
的 修正 提供 参考 . 


2 ”数据 和 方法 
2.1 ”数据 

美国 科罗拉多 大 学 基于 CHAMP 和 GRACE 了 
星 加 速 仪 资料 , 联合 GPS 接收 机 对 于 位 置 速 度 的 全 
弧度 测量 资料 进行 精密 定 轨 ,扣除 光 压 摄 动 后 ,1 
大 气 阻力 摄 动 模型 推导 出 热 层 中 性 气体 总 密度 , 数 
据 精 度 达 到 10- 攻 kg/ms 的 水 平 . 时 间 分 辩 率 为 458 
(CHAMP) 和 10s (GRACE). 本 文 还 采用 欧 空 局 
(ESA) 发 布 的 基于 SWARM-C 卫 星 加 速 仪 资 料 导出 
的 大 气 密度 数据 DNSxACC_2 (cdf 格 式 ， 网 址 为 


(高 度 、 地 方 时 、 太 阳 辐 射 等 ), 针对 相同 纬度 区 间 ， 
按照 30 准 则 计算 密度 比 的 平均 值 及 标准 偏差 , 最 终 
得 到 平均 密度 比 曲 线 , 可 识别 MDM 特征 , 提取 赤道 
附近 的 密度 局 部 峰值 和 中 纬度 附近 的 2 个 局 部 谷 值 . 

我 们 注意 到 , 日 下 点 的 位 置 变 化 会 引起 大 气 
密度 的 南北 半球 不 对 称 (夏季 半球 和 冬季 半球 的 差 
异 ), 对 识别 MDM 的 极 值 、 统 计 平 均 过 程 时 产生 外 
在 干扰 , 因此 尽量 选择 日 下 点 在 赤道 附近 区 域 的 资 
料 ( 太 阳 赤 纬 在 10° 以 内 ), 对 提取 MDM 特 征 更 有 利 ， 
而 当 资 料 不 充分 时 , 会 适度 扩大 到 15° 门 限 . 


3 ”结果 分 析 
3.1 “个 例 对 比 
为 展示 MDM 现 象 , 首先 选择 同一 天 的 两 颗 卫 


65 卷 


Ed LX 
星 资 


Z] 


料 做 对 比 ， 


不 同 地 方 时 的 观测 结果 . 第 1 组 
CHAMP 卫 星 落 在 夜间 半球 的 


BO Sita, 图 


1 给 出 两 颗 卫 星 轨 道 分 别 落 在 
2006 年 3 月 14 日 ， 
也 方 时 为 05:17 LT, 
1 (a) 没 有 明显 MDM 特 征 ; GRACE 


H 
KE 


ERKE: 高 层 大 气 密度 的 午夜 极 大 人 


特征 分 析 与 建 模 


2006 年 4 月 19 日 ，CHAMP 地 方 时 01 


前 后 ， 


方 时 为 20:40 LT, 图 


图 


1 (c) 有 清晰 的 MDM 现 象 ; 


VA S 


卫星 的 同时 测量 


1 (d) 未 呈现 MDM 特 征 . 对 上 
可 知 : MDM 现 象 只 发 生 在 


:54 LT, 在 午夜 
而 GRACE-A 地 


Li 


-A 地 方 时 在 23:11 LT, 属于 午夜 前 后 , 图 1 (b) 存 在 “午夜 前 后 (如 23:11 LT. 01:54 LT), 而 在 20:40 LT 和 
清晰 的 MDM 现 象 ( 赤 道 附近 的 峰值 ， 南 北半球 中 05:17 LT 则 没有 MDM 现 象 . 
纬度 地 区 的 两 个 谷 值 , 方形 阴影 区 所 示 ). 第 2 组 是 
-80 -60 -40 -20 0 20 40 60 -80 -60 -40 -20 0 20 40 60 80 
(a) date: 2006-03-14, 05:17 LT | 4 (b) date: 2006-03-14, 23:11 LT š 
l —a— CHAMP \ —A— GRACE-A À 
Ar T \ 414 
© 3r 7 A 3 
S 
azb % 十 412 
[7] 
$ “ck m T ail 
OF —m—w ^ T \ 11.0 
X a A ac 
L ee a M 并 4 
ý L L fi à L L fi L 1 V a 1 i 
4 (c) date: 2006-04-19, 01:54 LT (d) date: 2006-04-19, 20:40 A es 
3L —=— CHAMP 二 —^— GRACE-A |/^ 3 
A^ 
"d 
2r 十 1.2 
2 
© peak | 
E 1 " [ RS 4 1.1 
g A f ies m 
LA" EER 
valley 1 Li ^ 
0.9 - N PA 十 \ / 40.9 
A 
v NN / 
08r valley 2 + NS. 40.8 
-80 E X 20 0 20 40 60 -80 60 40 20 0 20 40 60 80 
latitude / ° latitude / ° 
图 1 CHAMP 和 GRACE-A 卫 星 在 不 同 地 方 时 的 同时 测量 及 其 密度 比 随地 理 纬度 的 变化 


Fig.1 The comparison of CHAMP and GRACE-A’s density ratio observation at different 


geographic latitude for same days 


3.2 
3.21 地 方 时 差异 


影响 MDM 特 征 的 因素 


采用 地 方 时 1 b 为 宽度 , 按照 22:00 一 23:00 LT. 
23:00 一 0:00 LT、0:00 一 01:00 LT、01:00 一 02:00 LT. 
02:00 一 03:00 LT、 03:00 一 04:00 LT 对 资料 进行 分 


ZH, 统计 得 到 不 同 纬度 


区 间 密 度 比 的 出 


上 线 ， 见 图 2. 


Ps 


每 个 紧 列 对 应 同一 颗 - 


结果 ， 自 上 而 下 按 地 方 


9-3 


时 区 间 分 成 6 行 , 取 各 


ocal time and its variation with 


自 中 值 22:30 LT. 23:30 LT. 


00:30 LT, 01:30 LT. 02:30 LT. 03:30 LT 做 标记 , 可 
以 看 出 : 在 22:30 LT (对 应 22:00 一 23:00 LT) 3 颗 卫 


0:00 LI) 开始 ，3 颗 卫星 都 记录 到 了 MDM 现 象 , 在 


Me 


AB 


星 都 没有 记录 到 MDM 现 象 , 总体 特征 是 高 纬度 区 
附近 最 低 . 从 23:30 LT (对 应 23:00 一 


区 域 出 现 密度 峰值 , 南北 半球 


纬度 出 现 谷 值 ， 


65 卷 


随后 的 0:00 一 01:00 LT. 01:00—02:00 LT 波 


向 南半球 移动 ; 
MDM 特 征 依然 显 


条 不 断 


在 02:00 一 03:00 LT，CHAMP 的 
著 , 而 GRACE-A 和 SWARM-C 记 


录 的 信号 开始 减弱 , 在 03:00 一 04:00 LT, MDM 完 全 


消失 . CHAMP 的 MD 


M 持 续 时 间 更 长 , 可 能 是 高 度 


差异 引起 的 ,间接 佐 
向 上 传播 所 致 加. 


证 MDM 可 能 是 低层 大 气 波 动 


为 定量 描述 MDM 特 征 的 强 弱 , 设计 “MDM 振 


Z] 


幅 * 统 计量 : 以 
同 纬度 的 密度 上 


2 (2-b) 为 例 对 于 2.2 
E, 求 -15° 一 15° 的 低 纬度 最 大 值 ( 记 


Hee AB 


"BART 


的 不 


Wp), 南半球 区 域 -25° 一 45° 最 小 值 ( 记 为 v1), d 


^E SR [X 325? —45^ zz [R] 
最 小 值 求 平均 , MAPA 
纬度 最 大 值 p 减 去 中 纬 


Br 
又 


9 一 2 


u 


的 最 小 值 ( 记 为 vz)， 将 两 个 


三 | 


度 密 度 的 


AR 
背景 值 , 得 出 MDM 振 幅 
Ui T V2 
2 


直 ， 再 用 低 


(1) 


比较 图 2 中 相同 卫星 不 同 地 方 时 的 振幅 , CHA- 
MP 和 GRACE-A 的 最 大 值 出 现在 0:00 一 01:00 LT, 
SWARM-C 的 最 大 值 出 现在 01:00 一 02:00 LT. 


CHAMP GRACE-A SWARM-C 
-80 -60 -40 -20 0 20 40 60 80 -80 -60 -40 -20 0 20 40 60 80 -80 -60 -40 -20 0 20 40 60 80 
1.20 L 1 T T T T T T 加 | 1.50 T T T T T T T T T T T T T 
135 [ E 
€ 1142F m J 135L A ^ 4 
a fae EN " a j 120 L N A J 1.20 上 oe、 .* 4 
3 "T N A 105 F s 7 4 
8 0.96 F TNR 4 105 L A a 4 doni e. 村 w 
0.88 È m (1-a) 4 0.90 È PEUT (2-a)4 0 元 让 ee (3-a) 
920 上 1 1 1 i 1 1 1 E 1.50 1 1 i 1 1 1 1 1 1 1 1 1 1 
"ow Po 870263, | 135E A 820.154 , 4 1.35 F 8=0.1137 4 
J 104} 、 JA P 120 上 上 入 ra 1.20 E e Z J 
9 oor EN / mo 4 1.05 L \ x 4 s F em 7 "E J 
~ osp — tV (1-5) 0.90 F Naa AQ AM (2-b)] 0.75 F “tee (3-b) 
oon 上 1 1 i fi 1 1 1 E 1.50 1 1 1 i L L . e 1 1 i 1 1 1 
1120 a /'w 8-0347. 4 135L ^ 8-0.250 4 135p X 8-0.149 e 7 
LE \ ron Y 人 120F v 4 
1.04 E x 4 .20 } J 
5 8 0.96 È K. / NN " 105 上 A A 人 A ^ J Hor ` ere 2 J 
w © 088 上 7 J a a X x“ 0.90 - e-e- Nex he J 
ki 0.88 LA (1-c) 0.90 - A^ a a^ (2-c)] 075 L ee (3) 
> 0.80 Et 1 1 1 1 1 LL 1.50 1 1 1 1 iA 1 1 1 1 1 1 1 
vaL Pe 5-020 | 1.35 F 8-0235 4 135E 9 5=0246 ] 
$ & P A 1.20F N e ] 
a 1.04 上 EN / " ad 120p A A 4 WEE. ^4 x X a 
8 0.96 F nu N "M or A AUR, ] oso E eu me Pa 4 
5 ossL i: Num (1-d) 4 0.90 È 从 A. A (2-d)4 ozsL Sue. (90) 
0.80 上 1 fi fi i fi fi fi 4 1.50 1 1 fi i LT 1 1 1 i «Leni L 
E [* | 4 
EROS / MNO 8=0.302 | 1.35 L J 1.36 F 8=0.128 。] 
I 14) a x p m] 120L 全 i n 120 e / 1 
3 oer Y / BN / 4 qos: 4*4 E 1.05 ^. * j 
a ° LE " "| : "a 090 L ooe Or J 
0.88 - Saa (9) 1 0.90 - "A-A. /2-e)4 ove ek "e^  (3-e)] 
0.80 - 1 1 1 i 1 1 1 4 1.50 1 1 fi i Aa L > L fi fi j 1 1 1 
1.20 F 4 ij 150F €. J 
5 1412F m 5=0.119 j 1.35 上 ^. 4 435L ee J 
g 104} Pa DS J 120b “a A J 1.20 E EM 4 
8 osf "m yoo d 1.05 È E A 4 108 F E^ PE 
0.88 È (1-f) 了 0.90 上 A iiA asa (2-f) 4 075 È LONE e-e-9(3-f) J 
0.80 E 1 1 1 fi 1 1 1 E 1 1 1 i 1 1 1 1 1 1 i 1 1 ji 
-80 -60 -40 -20 0 20 40 60 80 -80 -60 -40 -20 0 20 40 60 80 -80 -60 -40 -20 0 20 40 60 80 
latitude / ° latitude / ° latitude / ° 
图 2 3 颗 卫星 在 不 同 地 方 时 区 间 的 密度 比 随地 理 纬 度 变 化 的 统计 结果 (CHAMP, GRACE-A#ISWARM-C) 


Fig.2 The observation of density ratio with geographic latitude at different local time interval based on three satellites 
(CHAMP, GRACE-A and SWARM-C) 


3.22 ”高 度 差 异 


为 观察 不 同 高 度 大 气 的 MDM 有 无 差异 ， 


Z] 


3 给 


出 了 CHAMP 和 GRACE-A 在 同一 天 观测 的 个 例 


hi 


Xt Eb DA A IK EH 4 


较 . 首先 选择 2005 年 4 


射 情形 下 多 天 资料 的 统计 比 


10 


日 CHAMP 和 GRACE-A 
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的 资料 , 两 颗 星 地 方 时 相差 仅 为 0.01 h， 可 认为 几 
FHH, 高 度 在 374 km 和 481 km, 相差 110 km. 两 


颗 星 的 夜间 密度 曲线 如 
0.3586, GRACE- 


表明 高 度 越 高 振 


幅 越 小 . 


图 3 (a) 所 示 , CHAMP 振 幅 
A 振幅 0.3109, 后 者 是 前 者 的 87%， 
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图 3 不同 高 度 卫星 资料 的 密度 比 随地 理 纬度 变化 的 统计 结果 (CHAMP, GRACE-A) 


Fig.3 The observation of density ratio with geographic latitude at different altitudes based on CHAMP and GRACE-A 


步 选取 更 多 资料 做 统计 ， 资 料 筛 选 规则 
是 : 太阳 射电 辐射 流量 已 。; 小 于 100 sfu (1 sfu 
10-22 W-m-?-Hz~?), 地 方 时 23:00 一 02:00 LT, X 
阳 赤 纬 10" 以 内 . CHAMP 共 38 组 数据 ，GRACE-A 
有 103 组 数据 ， 图 3 中 (b) 和 (c) 同 时 给 出 平均 密度 比 
(点 线 ) 和 误差 棒 ( 浅 灰色 和 深 灰 色 阴 影 轮廓 )，CHL- 


AMP 振 幅 0.308，GRACE-A 振 幅 0.210， 是 前 者 的 
68%. 因此 从 统计 角度 同样 表明 高 度 越 高 , MDM 振 
幅 越 小 . 


表 1 给 出 3 个 大 气 模型 在 不 同 高 度 的 模拟 结果 ， 


反映 出 随 高 度 增加 振幅 减 小 的 趋势 , 但 在 具体 数值 
和 下 降 速率 上 与 观测 不 一 致 . 与 观测 对 比 结果 如 下 : 
(1) 在 360 km、480 km 两 个 高 度 上 , DTM2000 模 型 
的 振幅 仅 为 观测 的 46% 和 53%, NRLMSISE00 模 型 
仅 为 观测 的 33% 和 26%; (2) 两 个 高 度 之 间 的 差异 
(用 GRACE-A 所 在 高 度 的 振幅 除 以 CHAMP 所 在 
高 度 的 振幅 )， 观 测 为 68%、DTM 模 型 为 78%， 后 
者 低估 了 差异 ; NRLMSISE00 为 54%, 高 估 了 差异 ; 
(3) JB2008 则 完全 没有 反映 出 MDM 现 象 , 无 法 计 
算 振幅 . 


7 


表 1 不 同 高 度 观 测 和 3 个 大 气 模 型 模拟 的 MDM 振 幅 
Table 1 The MDM amplitudes of observation at different altitudes and simulation results by 3 


atmospheric density models 


Satellite Ave. Altitude Observation DTM2000 NRLMSISE00  JB2008 

CHAMP 360 km 0.308 0.143 0.101 / 
SWARM-C 450 km 0.170 0.126 0.057 / 
GRACE-A 480 km 0.210 0.112 0.055 / 


3.2.3 ”辐射 水 平 差异 


为 观察 不 同 太 阳 辐 射 水 平 下 MDM 有 无 差异 ， 
图 4 给 出 GRACE-A 在 辐射 差异 显著 的 两 天 观测 对 


比 以 及 3 种 辐射 水 平 的 统计 比较 . 首先 选择 两 天 
处 于 不 同 辐射 水 平 的 资料 : 2009 年 3 月 18 日 辐射 为 


9-5 


69 sfu ( 弱 辐 射 水 平 )，2012 年 9 月 25 日 辐射 137 sfu 
(中 等 辐射 水 平 ). GRACE-A 在 这 两 天 的 夜间 地 方 
时 很 接近 , 分 别 为 01:07 LT 和 01:18 LT. 图 4 (a) 显 示 
MDM 振 幅 分 别 为 0.594 ( 弱 辐 射 ) 和 0.101 (中 等 辐 
射 )， 后 者 仅 为 前 者 的 17%， 可 见 辐 射 水 平 越 强 
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MDM 振 幅 越 小 . 振幅 
选择 更 多 的 资料 进行 统计 研究 ， 围 绕 时 间 
跨度 最 长 的 GRACE-A 了 卫星， 选取 地 方 时 23:00 一 个 方面 的 原因 
02:00 LT、 太 阳 赤 纬 小 于 15° 为 总 样本 , 再 按照 辐射 对 振幅 的 影响 不 够 突出 


水 平 区 分 成 3 组 : 第 1 组 见 
第 2 组 见 


二 有 104 d% 
100 sfu， 共 有 45 d 资 料 ; 第 3 组 见 图 


zt El: 
E , 


于 100 sfu, 


结论 一 致 四. 


的 统计 结果 见 表 2, 同样 证 明 随 着 辐射 增强 ， 
弱 . 需 补充 说 明 的 是 : 表 2 中 SWARM-C (450 km) 的 


RE, MEWA 
0.127. 0.104), 


Lf 


> (1) 两 颗 星 高 度 差 异 30 kmJ 
; (2) 按 辐射 水 平分 组 之 后 ， 


比 GRACE-A (480km) 振 幅 小 , 与 上 一 节 * 高 度 
越 低 ，MDM 振幅 越 大 "的 结论 不 一 致 ， 可 能 有 两 


不 显著 ， 


图 4 (b), 辐射 小 于 80 sfu, 


19 d 资 料 , 从 图 4 (b)-(d) 的 统计 结 
寺 水 平 增 大 ， 振 幅 不 断 减 小 (0.219、 


E 


Hir 80- 情况 下 , 地 方 时 和 季节 等 其 他 因 


图 4 (c), 


4 (d), 辐射 大 


设 定 高 度 为 480 km, 


今后 需要 进一步 研究 . 
使 用 DTM2000、NRLMSISE00 和 JB2008 模 型 ， 
围绕 3 个 辐射 水 习 


与 Ma 等 


用 相同 分 组 对 CHAMP 和 SWARM-C 


基于 CHAMP 卫 星 资料 的 结果 见 图 


振幅 减 


5， 从 图 中 可 见 DTM2000 的 所 
0.117) 和 NRLMSISE00 振 幅 (0.0045-0.047) 随 fi 
水 平 变 化 不 显著 , 没有 反映 出 辐射 
趋势 ; JB2008 模 型 则 没有 刻画 出 MDM 现 象 . 


两 颗 星 的 样本 分 布 不 尽 相 同 , 在 高 度 差 异 不 明显 
素 干扰 了 统计 结果 ， 


的 


FE 进行 模拟 ， 
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Fig.4 The comparison of density ratio under different solar radiation levels based on the observation of GRACE-A 
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R fi (0.098— 


射 


越 强 振幅 越 小 的 
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图 5“ 大气 模 型 对 不 同 辐射 水 平 下 MDM 的 模拟 结果 (以 480 km 高 度 举例 ) 


Fig.5 The simulated results of MDM by three atmospheric density models under different solar radiation levels (take altitude 


480 km as example) 


表 2 对 于 不 同 辐射 水 平分 组 的 振幅 统计 结 
Table 2 The statistical results of MDM’s amplitude under different solar radiation levels 


CHAMP GRACE-A SWARM-C 
Solar radiation level 

Num. Amp. Num. Amp. Num. Amp. 

Fio.7 < 80 45 0.312 104 0.219 84 0.153 

80 < Fio.7 < 100 15 0.261 45 0.127 14 / 

Fio.7 2 100 27 0.182 19 0.104 21 0.059 
324 ”季节 差异 了 相应 的 模型 结果 , JB2008 模 型 未 刻画 出 MDM 现 
图 6 给 出 了 3 个 特殊 时 间 的 典型 个 例 对 比 : 第 1 Bs ERAT EA, E(b)MEE) D TM2000% 
行为 2007 年 6 月 6 日 , 临近 夏至 日 , 日 下 点 纬度 (图 中 ii 为 0.262、 0.284, NRLMSISE00 振 幅 为 0.113、 
标注 sub-solar lat.) 在 北纬 22.59"; 第 2 行为 2006 年 9 0-160, 最 着 大 于 秋分 日 (图 (qd)) 的 模型 振幅 ( 约 0.1)， 
用 22 日 秋分 日 , 日 下 点 纬度 在 赤道 附近 (北纬 0.46"); 。” 忆 CHAMP 观 测 所 示 的 “秋分 日 振幅 大 于 至 日 ”特征 
第 3 行为 2007 年 1 月 24 日 , 在 冬至 日 之 后 , 日 下 点 纬 ”恰恰 相反 . 因此 , 模型 对 季节 差异 刻画 存在 较 大 偏 

度 在 南 纬 19.359. CHAMP 卫 星 在 这 3d 的 阴影 区 地 ” 差 . 

方 时 分 别 为 23:54 LT、23:35 LT 和 0:08 LT, 相互 之 选择 更 多 天 的 资料 ,按照 太阳 赤 纬 的 范围 进 
间 很 接近 . TAL, 辐射 流量 互 0.7 指 数 都 在 80 sfu 左 行 分 组 统计 , 见 表 3. 第 1 组 是 太阳 赤 纬 在 北纬 15°- 
右 ,地 磁极 其 平静 (日 均 地 磁 指 数 Ap = 2), 轨道 23° (夏至 日 前 后 ); 第 2 组 是 太阳 赤 纬 在 10° 以 内 
高 度 在 358-368 km 之 间 ， 图 6 左 列 是 CHAMP 观 测 (EA KOME); 第 3 组 是 太阳 赤 纬 在 南 纬 15°- 
得 到 的 密度 比值 , 图 (a)、(c) 和 (e) 都 展示 出 清晰 的 ”23”( 冬 至 日 前 后 ). 选取 地 方 时 23:00 一 02:00 LTH. 
MDM 现 象 , 峰值 位 置 始终 位 于 夏季 半球 低 纬度 (图 0.7 小 于 100 sfu 的 资料 , 样本 最 充分 ， J 
(a) 的 北纬 5°, 图 (e) 的 南 纬 15°), 比较 秋分 日 (图 (c)) 至 点 附近 (夏至 、 冬 至 ) 振 幅 普 遍 比 分 点 (春分 、 秋 
和 至 日 (图 (a)、(e)) 的 振幅 ,发现 秋分 日 的 振幅 大 分 ) 小 , 仅 为 后 者 的 42%-52% (CHAMP), nene 
于 两 个 至 日 的 振幅 (2.6 倍 和 2 倍 )， 图 6 的 右 列 给 出 (GRACE-A). 2596-7296 (SWARM-C), 说 明 季 节 
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效应 削减 了 MDM 引 起 的 振幅 ; 总 体 看 ， 峰 值 位 置 GRACE-A 在 夏至 、SWARM-C 在 秋分 时 , 峰值! 


上 上 


大 多 出 现在 南半球 (-5 和 -15")， 只 有 CHAMP 和 现在 北纬 5 附近 . 


Fig.6 The 


Table 3 
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图 6 季节 效应 对 MDM 的 影响 : 夏至 、 秋 分 和 冬至 附近 3 d 的 CHAMP 观 测 和 模拟 结果 
seasonal effect on MDM: the CHAMP’s observation and models’ simulated results on three days around summer 


solstice, autumn equinox and winter solstice 


R3 不 同 太 阳 赤 纬 时 MDM 振 幅 和 峰值 位 置 的 统计 结果 
The statistical results of MDM’s amplitude and peak position with different solar latitude 


regions 
CHAMP GRACE-A SWARM-C 


Solar Latitude 
Amp. Peak pos. Amp. Peak pos. Amp. Peak Pos. 


15°—23° 0.128 5° 0.101 5° 0.042 —5° 
—10°—10° 0.308 —5° 0.210 —5° 0.170 5° 
—23?——15? 0.161 —5° 0.124 —5° 0.122 —15° 
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4 MDM 特 征 经 验 建 模 

如 前 文 所 述 , 大气 模 型 对 MDM 的 刻画 存在 明 
EDE, 本文 尝试 基于 3 颗 卫 星 23:00 一 02:00 LT 的 
MDM 信 号 , 建立 有 效 的 经 验 函 数 形式 对 观测 进行 
WA, 力求 更 准确 地 刻画 MDM 特 征 . 经 反复 试验 ， 
先 选择 地 理 纬度 的 6 阶 勒 让 德 函 数 作为 基础 ， 


yle) = ao + oP(p) + a Pa(p) + az3Ps(u)+ 
a4 P4(u) + asPs(u) + asPs(u). 

再 加 入 与 地 方 时 、 高 度 的 耦合 项 

y(p,t, Z) = ao + [a1 Pı (u) + as P»(p) + as Pau) 9- 
a4P4(p) + a5 Ps + ag Ps(u)|+ 
t- [by Py (u) + b2 Po(w) + b3 Ps(u)+ 
baPa(u) + bs P5 (tu) + bs Ps(u)]9- 
Z le Piu) + e P2() + e3 Py(p)9- 
ca4Pa(u) + cs P5(u) + csPs(u)], — (3) 


其 中 4 = sing, 2 是 地 理 纬 度 , P(1) 是 第 i 阶 勒 让 德 
函数 , 本 文 取 ; = 1,2,--- ,6, t 是 地 方 时 , Z 是 高 度 ， 
ao, i, bi, Ci 是 拟 合 系数 . 
图 7 给 出 了 用 不 同 经 验 函 数 对 MDM 观 测 特征 
的 拟 合 情 况 . 对 于 单一 纬度 因素 的 函数 (2) 式 ， 拟 
合 结果 如 图 7 (a)， 相 关系 数 "为 0.738, 标准 偏差 为 
0.08, 展示 了 清晰 的 MDM 特 征 : 峰值 在 南 纬 5°， 两 
个 谷 值 在 南 纬 45*、 北 纬 35*, 整体 偏向 南半球 ， 中 
纬度 谷 值 0.91， 赤道 附近 峰值 1.17, 振幅 为 0.26, 相 
当 于 在 平均 密度 为 1 的 比例 尺 下 , MDM 引 起 中 低 纬 
度 密度 大 约 26% 的 起 伏 . 

考虑 地 理 纬 度 因素 与 高 度 、 辐 射 、 地 方 时 的 
耦合 效应 ((3) 式 ) 后 ， 整 体 拟 合 效果 得 以 显著 改善 ， 
7 为 0.923, 标准 偏差 为 0.046. 分 别 对 每 颗 星 画图 和 
统计 ， 图 7 (b) 显 示 CHAMP 资 料 拟 合 得 最 好 , 7 为 
0.972; 图 7 的 (ce) 和 (d) 分 别 是 GRACE-A 和 SWARM- 
C 的 拟 合 结 果 , 7 分 别 为 0.868 和 0.886. 横 轴 是 样本 序 


(2) 
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号 , 分 成 3 组 地 方 时 , 拟 合 曲线 与 观测 符合 得 较 好 ， 
且 刻 画 了 不 同 峰值 谷 值 位 置 以 及 与 地 方 时 、 高 度 
等 的 关系 . 


5 结论 

利用 CHAMP、GRACE-A 和 SWARM-C 等 3 颗 
极 轨 卫 星 的 资料 , 研究 了 大 气 密度 在 低 纬 度 区 域 的 
午夜 极 大 值 (MDM) 现 象 . MDM 一 般 出 现在 地 方 时 
23:00—02:00 LT 之 间 , 峰值 位 置 在 低 纬 15" 以内, 谷 
值 位 置 在 中 纬度 35°-45°, 整体 更 偏向 南半球 , 振幅 
为 平均 密度 的 26% 左 右 . MDM 同 时 受到 高 度 、 辐 
射 水 平和 季节 的 影响 , 随 着 高 度 增 大 及 辐射 水 平 的 
增强 , MDM 振 幅 均 呈 减 小 趋势 ， 冬至 夏至 前 后 的 
MDM 振 幅 较 弱 ， 春 秋分 附近 的 MDM 振 幅 最 大 . 
将 这 些 观 测 特征 与 DTM2000、NRLMSISE00 和 JB 
2008 模 型 的 模拟 结果 对 比 , 发 现 : (1) JB2008 没 有 
刻画 出 MDM 现 象 ，(2) 模 型 低估 了 MDM 振 幅 ， 在 
360 km, 480 km 两 个 高 度 ,，DTM2000 模 型 的 振幅 
4X 2331130 (5146968815396, NRLMSISE00 模 型 仅 为 观 
38339612696; (3) 模 型 没有 准确 刻画 MDM 与 高 
度 、 辐 射 水 平和 季节 的 关系 . 
基于 3 颗 卫 星 观测 资料 联合 求解 ,研究 了 一 种 
基于 地 理 纬 度 的 6 阶 勒 让 德 多 项 式 ， 同 时 融合 地 方 
时 和 高 度 因素 的 经 验 函 数 , 较 好 地 描述 了 MDM 特 
TE, 相关 系数 为 0.923, 标准 偏差 0.046, 显著 优 于 仪 
考虑 纬度 因素 的 函数 形式 (相关 系数 0.738, 标准 偏 
250.080), 拟 合 曲线 在 振幅 和 相位 上 与 观测 符合 得 
BUT. 该 经 验 函 数 可 为 大 气 密度 模型 修正 提供 参考 ， 
致谢 感谢 美国 科罗拉多 大 学 提供 了 基于 CHAMP 
和 GRACE 卫星 加 速 仪 资料 反 演 的 大 气 密 度数 据 
(https:/ /lasp.colorado.edu/space-weather-portal); 
感谢 欧 空 局 ESA 提 供 SWARM 了 卫星 的 数据 (https: 
//earth.esa.int /eogateway /missions/swarm/data, 


ftp:/ /swarm-diss.eo.esa.int). 
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Fig.7 The comparison of MDM observation and their fitting curves based on two kinds of empirical function 
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The Feature Analysis and Modeling of Upper Atmospheric 
Midnight Density Maximum 
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(2 Key Laboratory of Space Object and Debris Observation, Chinese Academy of Sciences, Nanjing 210023) 


Asstract The features of upper atmospheric midnight density maximum (MDM) around low geo- 
graphic latitudes are studied based on neutral mass densities data at altitudes 360—480 km, derived from 
the accelerometer measurements aboard on the three polar orbiting satellites CHAMP (CHAllenging 
Minisatellite Payload), GRACE-A (Gravity Recovery and Climate Experiment-A) and SWARM-C (The 
Earth’s Magnetic Field and Environment Explorers-C). The MDM appears during the local times from 
23:00 to 02:00 LT, whose peak locates at the low latitudes within 15° and two valleys locate at the mid- 
dle latitudes between 35° and 45° on both hemispheres separately. The structure of MDM drifts toward 
the southern hemisphere overall. The MDM's amplitude decreases with increases in altitude and solar 
radiation level. The seasonal effect weakens the MDM's amplitudes around the summer and winter sol- 
stices, while the amplitudes around the spring and autumn equinoxes are extremely significant due to the 
slight seasonal difference between both hemispheres. Three atmospheric density models DTM2000 (Drag 
Temperature Model 2000), NRLMSISEO00 (Naval Research Laboratory Mass Spectrometer and Incoherent 
Scatter Radar Extended atmosphere model) and JB2008 (Jacchia-Bowman 2008 model) are used to sim- 
ulate the MDM along these three satellites’ orbits, and compared with the observations. It is found that 
the JB2008 model is failed to describe the MDM, and the other two models underestimate the MDM’s 
amplitudes at altitudes 360 km and 480 km: the simulated amplitudes by the DTM2000 model are 46% 
and 53% of the observed amplitudes respectively, and only 33% and 26% for the NRLMSISEO00 model. 
These three models are also failed to depict the MDM’s variation with altitude, solar radiation level and 
seasonal effects. In order to correct the model prediction, a 6th-order Legendre polynomial of geographic 
latitude, coupled with arguments of local time and altitude, is designed to fit the MDM signals from 
the three satellites’ observations. In terms of amplitude and phase of the MDM, the fitting results agree 
with the observations very well, and the correlation coefficient is 0.923. It indicates that this empirical 
polynomial could be helpful to the density model correction and high accuracy prediction of spacecrafts 
in low Earth orbits. 


Key words atmospheric effects: atmosphere density model, methods: data analysis 


